Universal Properties of Ferroelectric Domains 
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Basing on Ginzburg-Landau approach we generahze the Kittel theory and derive the interpolation 
formula for the temperature evolution of a multi-domain polarization profile P(x, z). We resolve the 
long-standing problem of the near-surface polarization behavior in ferroelectric domains and demon- 
strate the polarization vanishing instead of usually assumed fractal domain branching. We propose 
an effective scaling approach to compare the properties of different domain- containing ferroelectric 
plates and films. 

PACS numbers: 77.80.Bh, 77.55.+f, 77.80.Dj 



Design of ferroelectric devices necessitates taking into 
account such finite size effects as the formation of 
polarization-induced surface charges that, in turn, pro- 
duce the energy consuming electrostatic depolarizing 
fields (see Ref.[l] for review). As a result, regular pe- 
riodic structures of 180^ domains that alternates the sur- 
face charge distribution, firstly proposed by Landau and 
Lifshitzjiflsj and by Kittel [1| for ferromagnetic systems, 
can be formed in uniaxial easy-axis (natural or stress- 
induced) ferroelectric plates or films as an effective mech- 
anism to confine the depolarization field to the near- 
surface layer and reduce its energy (FigH^). The energy 
balance between the field-penetration depth (^ domain 
width d) and domain wall (DW) concentration (~ d~^) 
leads to the famous square-root Kittel dependence of d 
on the film thickness 2a f ^MMM- 



21C(3)-^-^^' W 

where ey and e±_ are the longitudinal and transversal 
dielectric constants and ^oa? is the transverse coherence 
length (roughly equal to the DW thickness). 

Consider the standard geometry [5] when the uniaxial 
ferroelectric film is sandwiched by electroded paraelec- 
tric passive layers of width and permittivity £p. The 
multi-domain state should exist in certain intervals of 
film thickness 2a f as shown in phase diagram in Fig {It 
and defined by the condition that delineates the applica- 
bility of Eq.([T]) and of our further consideration: 



^Ox < d(2af) < ap 



(2) 



the dependence d{2af) being given by ([T]). We also as- 
sume the most realistic case <C ^ ^ii ^^^^ gives 
d <^ af. At this stage the properties of domain struc- 
ture do not depend on a^, Sp and electrodes. For thicker 
films, when d{2af) approaches to ap the emergent depo- 
larizing field interacts with screening electrodes, Eq.ffl 
is not valid anymore, d growth exponentially with [5| 
and domains practically emerge from the sample. How- 
ever in free standing electrodless sample {ap oo) Kittel 



domains can exist in a wider interval of 2a f unless an- 
other restricting mechanism of the internal free charges 
screening does not came into the play. For thinner films 
we are turning to the region of little-studied atomic-size 
(microscopic) domains [8]. 

While domain structures should play a crucial role in 
the properties of thin ferroelectric films, only a few theo- 
retical analytical studies of their temperature dependence 
have been performed. In particularly the mostly used 
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FIG. 1: (a) Multi-domain texture of ferroelectric polarization 
in uniaxial ferroelectric film, sandwiched by two paraelectric 
(dead) -layers. The emerging depolarization electric field is 
provided by alternating polarization-induced surface charges 
and confined in the near-surface layer of thickness, compara- 
ble with domain width d. (b) Elliptical functions y = sn(x, m) 
for different parameters m that we use to model the domain 
profile at different t. (c) Phase diagram of domain states as 
function of sample thickness 2a/ and reduced critical tem- 
perature t — T/TcQ — 1. Polarization profiles of hard and 
soft domains were obtained by numerical solution of equations 
dSJ-dH). We assume that k\\ 500, £± 100, Sp <^ £±,>^||, 
^oa; — Inm and ap ^ 30nm. 
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Kittel approach 0, 0, in which the domain texture is 
considered as a set of up- and down- oriented (hard) do- 
mains, having a flat polarization profile P{x^z) = iPo^ 
DW are supposed to be infinitely thin and boundary 
effects on the ferroelectric-paraelectric interface are ne- 
glected, is valid only far below the transition temperature 
Tc. Although the more general consideration, proposed 
by Chensky and Tarasenko (CT)[9] (see also [S, J^]) is 
based on Ginzburg-Landau equations coupled with elec- 
trostatic equations is valid in the whole temperature in- 
terval, only the solution close to Tc was found. 

It is the objective of the present communication to es- 
tablish the approach that permits to model the tempera- 
ture evolution of domain structure. Basing on CT equa- 
tions we derive the analytical expression ([19]) for domain 
polarization profile that is valid in the whole tempera- 
ture interval and includes the Kittel (at low T) and CT 
(at T = Tc) solutions as particular cases. Then, we de- 
duce universal scaling relations between parameters of 
the multi-domain state that should be useful in treatment 
of experimental data. Our approach is complimentary to 
the frequently used first-principia simulations (see e.g. 
fllj). that reproduce the domain structure but give no 
general vision and parameter dependence of the results. 

Deducing the CT equations we are basing on the Euler- 
Lagrange variational formalism, that permits also to ob- 
tain the correct boundary conditions as variation of sur- 
face terms. The generating energy functional is written 
as S: 



taking into account the depolarizing field: 



/• 



F = I $(P, F,)dxdz, $(P, E) = $(P, 0) - EP - -^E^ 

OTT 

(3) 

where E = [E^.E^), P = [Px.Pz) and the field- 
independent part 

$(P,0) = ^>| + ^i/^. + ^/(P) (4) 

£±_ 2 £i\\ 2 K\\ 

includes the transversal P^, and non-polar longitudinal 
Pzi noncritical contributions {£^^£i\\ ^ 1). The nonlinear 
Ginzburg-Landau energy depends on the spontaneous z- 
oriented polarization P (assuming that = Pzi + P) 
and is written as: 



/(p) = + ip-2p4 - 



%aPf + ^(5.P)' (5) 



where the reduced temperature t is expressed via the bulk 
critical temperature as: t = T/Tcq — 1, parameter k\\ 
is expressed via paraelectric Curie constant C and via 
longitudinal zero-temperature permittivity £\\ in ([T]) as: 
X|| = C /TcQ — 26: II , and coefficient Pq is roughly equal 
to the saturated bulk polarization at T <C Tc 

The variation of (|3]) with respect to P and the elec- 
trostatic potential (E = —Vip) and excluding of the 
non-essential variables Px and Pzi gives the system of re- 
quired equations that describe the ferroelectric transition 



- ^Ldl - e,.dl)P + (P/Po)^P = -^dz^. (6) 

{£,\\dl^£^dl)^ = A^dzP 

These equations should be completed by the Poisson 
equation for paraelectric media in which ferroelectric film 
is embedded: 

{dl+dl)ip^^^=Q, (7) 
and by boundary conditions at the Para-Ferro interface 

£i\\dz^ - £pdz^^''^ = 4^P, ^ = (^^^\ dzP = 0. (8) 

that are also obtained as result of variation of (|3]) [l5[. 
Periodic conditions 

P(x, z) = P{x + 2d, z) (p{x, z) = (p{x + 2d, z) (9) 

with variational parameter d are imposed to describe the 
periodicity of domain structure. 

A simplification can be achieved if present the initial 
functional (|3]) using the dimensionless (prime) variables: 

z = afz\ X = r"^/^^oa^x , t = rt\ (10) 
1 



P = rl/'PoP^ (p= — r^/'a/Po(/9^ 



r^/^Po^P' 



with 



r=(^V^«l 



(11) 



in truncated form. 



F' = y■[4^(^i^'P'2 + ip'4+l(5;p'^2 



Stt 



{d'.^'f + P'd',^']dx'dz' 



(12) 



that was obtained after neglecting the small terms 



(13) 

(justification is given in Appendix) and minimizing over 
P P 

The Euler-Lagrange variation of ^T2\i over P' and cp^ 
gives the corresponding dimensionless equations: 



(14) 

5^^' = 47ra;P', (15) 
and boundary conditions at z'_ = and at z'_^_ = 2a^ = 2: 

P' = 0, = (16) 
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that are simpler then conditions ([8]) since the order of (|6]) 
was reduced by neglecting (p!3|) . We stress here that these 
conditions are derived from functional ([12]) as variational 
surface terms. 

Passage to dimensionless variables is the powerful tool 
that permits to study the various properties of ferro- 
electric domains even without solution the differential 
equations. Note first that equations (|14|15p contain only 
one driving variable - the dimensionless temperature f . 
Therefore the "master" temperature dependence of any 
physical parameter calculated from (|14|15p can be re- 
scaled for any other ferroelectric sample, using the rela- 
tions (fTQ|) . 

We derive now such "master" variational solution of 
equations (|14|15p for domain profile P' {x' ^ z' ^t') valid in 
the whole temperature interval. Note first that these 
equations can be solved analytically close to the tran- 
sition to a multi-domain ferroelectric state [1, [13] that 
occurs at: 



t' = — TT, tn = — 27r 



^11 ^Qx 

£^ 2a f 



(17) 



(in dimensionless and dimensional variables), when po- 
larization has the sinusoidal (soft) distribution: 



(18) 



with the half-period d'^ = \ph^ (that is expressed as ([T]) 
in dimensional variables but with 7 = tt and ey = X||/2). 
At lower temperatures domain walls become sharper due 
to the admixture of higher harmonics. At lower temper- 
atures the domains recover the (hard) Kittel-like profile. 

To account for both these cases by the unique inter- 
polation formula we shall exploit the depicted in Fig. [TJd 
periodical elliptical sinus function y = sn(x, m) = sn(x + 
4K, m) , frequently used to describe the incommensurate 
phases [12|. The 1/4 of the elliptical sinus period is given 
by the tabled first kind elliptical integral K{m) fl3]. The 
useful property of sn(x, m) is that, depending on the pa- 
rameter < m < 1 it recovers the all described above 
domain regimes: from the soft one (p!8|) at m = when 
sn(x, m) sinx (like in Eq. [T8|) to the hard (Kittel-like) 
one at m ~ 1 when sn(x, m) step-wise function. 

After some algebra (justification is given in Appendix) 
we arrive to the following variational expression: 



P' = A{t') sn 



2d'(t') 



x' ^ mi{t') 



sn [K2{t')z',m2{t') 



(19) 

where the temperature dependencies of parameters m\ it) 
and m2(t), elliptic integrals Kiit) and i^2(^), amplitude 
A{t) and domain lattice half-period d{t) are presented in 
Fig|2] and for practical use are approximated as: 



A\t') tanh0.35(t'-tO, 



d'{t') c^2.6 (20) 
tanh0.27(t'^ - t') 
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FIG. 2: Temperature dependencies of parameters of Eq.([T9|): 
(a) elliptic arguments mi and 7712, elliptic integrals Ki and 
K2 , (b) domain amplitude A and domain lattice period d' . 
All the variables are dimensionless. 



Formula (p!9|) satisfies the boundary conditions 

P'{x', z') = P'{x' + 2d', z'), P\x\ 0) = P\x\ 2) = 0, 

(21) 

recovers the soft domain structure ([18]) at t'^ when 
mi2{t'J = 0, A{f) - {f^ - ffl^ and the Kittel-like 
structure at low t'^ when mi2{t'^) 1, A{t') {—t'Y^'^, 
and gives the domain profile at arbitrary t' . Parameters 
Ki2{t') determine the space scale of polarization varia- 
tion: in dimensional variables the characteristic domain 
wall thickness is S,x{t) = i^x/ {—'^Y^'^ whereas the thick- 
ness of the near-surface layer where P{z) restores its equi- 
librium value is ~ d/{—tY^'^ • {>c\\/e\_Y^'^ (i.e. ~ d at low 

t). 

Variation and vanishing of polarization at the sam- 
ple surface modifies the initial assumption of the Kittel 
model that polarization is permanent inside domain and 
resolves the long-standing paradox [13] according to 
which the permanent domain polarization should be re- 
oriented close to sample surface by its own depolarization 
field that exists in the near-surface layer. 

As it follows from our calculations, the nonuniform dis- 
tribution of polarization pumps the depolarization charge 
p{r) ~ divP from the sample surface inside the near- 
surface layer ~ d, reducing the unfavorable depolariza- 
tion field (justification is given in Appendix) and its en- 
ergy £d ^ 1^.11 ~ 47rP^. The price of this - the 
dumping of the condensation energy Ec ^ 4.'KP'^/e\\ is 
not so high because ey ^ 1. That's why we believe 
that the near-surface polarization vanishing is more ef- 
fective mechanism to overcome the Kittel paradox in fer- 
roelectrics and reduce the near-surface depolarization en- 
ergy then the usually assumed [1, [13] but rarely observed 
fractal branching of alternatively oriented permanent- 
polarization domains near the sample surface. 

Polarization decay at the surface is the consequence 
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a) ^ b) ^ c) ^'^^^'^ 

FIG. 3: Polarization of Kittel domain, (a) Numerical solution 
of complete CT equations dn])-®. (b) Numerical solution of 
simplified equations (|14|) - (|16p . (c) Interpolation formula (|19p 

of the boundary condition = of simplified equations 
([HI)- ([16]). The validity of this effect is illustrated in Fig. [3] 
where we compare the numerical solution of simplified 
equations ([n|) -([T6 |) (Fig. [8)3) with that for the complete 
set of CT equations (Fig. |3^). Clearly the tendency of 
polarization vanishing is conserved for the case of gen- 
eral solution in Fig. [3^, although the "real" boundary 
condition dzP = ([8]) is satisfied exactly at the surface. 
Interesting to note that the precursor of the competi- 
tive surface domain branching is also seen at Fig. |3^b 
as ripples at the domain end-points. The corresponding 
variational solution ([T9|) at Fig. |3t is more smooth, but 
correctly represents the properties of numerical profile. 

We present now several remarkable conclusions about 
the physical properties of the multi-domain state which 
can be obtained only from the scaling properties ([TQ|) , 
without solution of CT equations ([6])- ([9]). 

(i) Any transverse length parameter scales as r~^/^^oa?- 
This, in particular, justifies the Kittel formula ([1]) for the 
domain width d even beyond the flat domain approxima- 
tion. A convincing demonstration of the validity of this 
scaling law was reported recently for various ferroelectric 
and ferromagnetic materials [7]. The temperature de- 
pendence d{t) can be incorporated into ([1]) as dependence 
7 = j{t). Meanwhile, the results shown in Fig. [2)3 as well 
as finite-element simulations [6] indicate that the depen- 
dence d{t) is very weak and hence one can extend the 
parameter 7 :^ 3.53 from ([T]) to any temperature. This, in 
particular, implies the low temperature hysteresis related 
with motion of DW. 

(ii) The temperature t scales as r. Thus, to com- 
pare the domain-provided physical properties of different 
plates or films (even constructed from different materi- 
als) it should be instructive to trace their temperature 
dependencies using the re-scaled coordinate t/r. 

(iii) All the domain-related properties and, in partic- 
ular, the transition temperature tc (flTj) and the soft- 
to-hard domain crossover temperature t* ~ lOtc scale 
as l/2a/ with plate (film) width, as illustrated in 
Fig. (it. The temperature interval for the existence of 
soft-domains At = tc — growth dramatically with de- 
creasing film thickness and one can expect that for thin 
films with 2a/ < lOOnm only soft domains with a gradual 
polarization distribution are possible. 



Summarizing we conclude that domains in any ferro- 
electric sample and at any temperature can be easily ob- 
tained from interpolation formulas (|19|2Qp applying the 
scaling relations ([TQ|) . This can be especially helpful to 
treat the experimental data, involving the local field dis- 
tribution of polarization inside domains like ESR or Ra- 
man spectroscopy, TEM domain imagery etc. 

We demonstrated that depending on the temperature 
and sample width domains can have soft (gradual) or 
hard (Kittel) profile. In any case polarization has the 
tendency to vanish at sample surface. 

Basing on universal scaling relations ([TQ|) we have 
demonstrated how the physical properties of the differ- 
ent multi-domain films can be compared and mapped 
onto each other. We hope that such method will give the 
power tool for analysis and systematization of numerous 
experimental data for thin ferroelectric films. 

This work was supported by the Region of Pi- 
cardy, France, by STREP "Multiceral" (NMP3-CT-2006- 
032616) and by FP7 IRSES program "Robocon". We 
thank to Prof. M. G. Karkut for the useful discussions. 



APPENDIX (EPAPS document) 

We present here the technical derivation of (i) sim- 
plified equations and corresponding boundary conditions 
from the generating Euler-Lagrange functional (ii) inter- 
polation formula for domain polarization (iii) justifica- 
tion of simplification of generating functional. 

We use the defined in the article dimensionless vari- 
ables, omitting the prime index. 

(i) Derivation of simplified equations and boundary 
conditions from the Euler-Lagrange functional 

Euler-Lagrange variation of the simplified dimension- 
less functional ([T2|) that describes ferroelectric phase in 
a infinite thin plate (film) located at < z < 2: over 
polarization P and potential of electric field cp gives: 

^ / r 4^ {tP5P + p^5P + (s,p) (5,jp) 1 

J [ -i^{da:^){dj^) + SPd, + SPd,^ + Pd,S^ \ ' 
An JSP {tP + P3 - dlP + -^d^i^) dxdz 
- -j5^{i^dl^-d,P)dxdz^[JP6^dx]lZl 

Two first (volume) terms provide the corresponding di- 
mensionless equations (fT4|) - ([T5]) whereas the third (sur- 
face) term gives the boundary condition ([T6|) that should 
be completed by condition of continuity of potential at 
z = and at z = 2. 
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(ii) Derivation of interpolation formula for domain po- 
larization 

Although the nonhnear equations (p!4|) - (p!6|) can not be 
solved exactly we shall look for their x-periodic domain 
solution in the variational form 



P 

/(O) 



/(^)sn 



4i^(mi) 
2d~~ 



■ X, mi 



(22) 



/(2) = 0, P{x,z) = P{x^d,z) 



considering mi, d and the function f{z) as variational 
parameters that minimize ([12]) 

Substitution of ([22]) back into ([12]) and integration over 
domain period gives: 



--^{dx(p)'^ + Pdz(p]dxdz 

OTT 

113 1 



(23) 



27r 



, ^^{da:(pf - ^ipdlip]dxdz 
J {djf {2dfS{rm)dz 



and 



47r I itp2 + 



4-K 



it/(^)2a(mi) + i/(^)4ry(mi) 



Now the functional depends only on variable z: 



47r 



1 (a(mi)t+ ^§^p{m,)) f{z)^ 



dz 



(24) 



where the coefficients are expressed via complete elliptic 
integrals of the first and second kind K{m), E{m) as: 



aim) 
r]{m) 

S{m) 



< sn^(x, m) > 
E{m)' 



1 

m 



1 



/3(m) 



K{m) 
< sn^(x, m) > 

— [2(1 + m)a(m) - 1] 
3m 

(S^ (x,m)) 
[4i^(m)]^ 
8 

m [4i^(m)] 
4:K{m) ((sn'ix) 



(25) 
(26) 

(27) 



E 

1=1,3,5 
2 



1 q^^'^{m) 



11 — q\m) 



4.K{m)- [2 - (1 + m) a{m)] 

o 



(28) 



Here q{m) = e ^("^) ^, S(x,m) = / sn{u,m) du and 
(. . .) is the average over the period. 
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FIG. 4: Coefficients a(m), /3(m), 7(m) and 5(m) that enter 
into the variational functional ()24p 



Dependencies a(m), /3(m), 7^(m) and (5(m) are pre- 
sented in FiglU 

The variational Euler-Lagrange minimum of ([M]) is 
given by the function: 



/(z) = A{t, mi, m2) sn [i^(m2) z, 1712] (29) 



with 



A(t,mi,m2) = 2d{t,mi,m2)K{m2) W2^^^m2 (30) 

that matches the boundary conditions /(O) = /(2) = 
providing that the dependence (i(t,mi,m2) is fixed by 
biquadratic equation: 

{2d)^S{mi){l^m2)K\m2)^{2d)^a{mi)t _ 

+4i^(mi)/3(mi) ~^ ^'^^^ 

Substitution of ([29]) back into ([24]) gives: 

F(mi,m2) = -47r^r]{m,) J f\z)dz (32) 

= -^7r7/(mi)A'^(t,mi,m2)7?(m2) 

Collecting all the results, we present the final variational 
solution ([T9D. 



(Hi) Justification of simplification of generating func- 
tional 

The simplified functional ([T2]) was obtained by neglect- 
ing the terms ([T3]) . 

Using profile P{x^z) from ([T9]) we can now justify that 
contribution of these terms is indeed small by noting that 
their action is concentrated in the near-surface layer of 
thickness ~ l/i^2(^) ^ We will consider only 

the Kittel regime far from tc = — tt. The soft regime close 
to tc was already considered in [16] . 

The relative contribution of the first term to F is esti- 
mated as 

/ Mdxdz/F ~ (£i)i/2^2i^^ ^ {- f\t\-'^' « 1 

(33) 
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that is small for the Kittel domains with d <C aj. 
Note, however, that this criteria is not satisfied for mon- 
odomain polarization profile that formally is achieved 
when d ^ oo. This means that dimensionless equa- 
tions (|14|15p can not be applied for monodomain x- 
independent solution, that however is unstable towards 
domain formation anyway. 

Another term A2 is related with the energy of the de- 
polarizing electric field Ez- According to (p!5|) . this field 
can be calculated from the polarization profile (p!9|) as: 

Ez{x^z) = —dz(p = — 47r9^ / / P{x2j z)dx2dxi. 

(34) 

It follows that the depolarization field Ez periodically 
alternates in x— direction in anti-phase with P and is lo- 
cated in the near-surface layer of thickness ^r- It vanishes 
at the surface and in the bulk. Estimating the maximal 
value of Ez at x ~ as ^^^^nax ^ we have: 

/ A2dxdz/F ~ fl(!l)l/2i2^^^ ^ f^(A)2|i|-l/2 « 1. 

(35) 

The physical meaning of this estimation is discussed in 
the main text of the article. 
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